library(ggplot2)
#dyn.load("test/Rprobit.so")
#dyn.load("test/IBIS.so")
#dyn.load("test/Rprob_separated.so")
dyn.load("test/Rinterface.so")
#dyn.load("test/ModelSelection.so")
#dyn.load("test/RLogit.so")
#dyn.load("test/Rmoral.so")
#dyn.load("test/RIBIS.so")
source("~/Documents/Memoire/src/Generate.R")
#source("generate_prostate.R")

M<-10000
ft<-function()
{
#p<-8
#n<-200
#library(MASS)
#X<-as.matrix(Pima.tr[,1:7])
#Y<-as.numeric(Pima.tr[,8])
#Y[Y==1]<-0
#Y[Y==2]<-1
#X<-X
#v<-apply(X,2,sd)
#m<-apply(X,2,mean)
#X<-t(apply(t(apply(X,1,"-",m)),1,"/",v))
#X<-cbind(1,X)
#source("R/Generate_Credit.R")
p<-ncol(X)
n<-nrow(X)
theta<-matrix(0, nrow=M,ncol=p)
w<-matrix(0, nrow=M,ncol=1)
t=proc.time()
res1<-.C("Maruyamaf", as.integer(M),as.integer(n),as.integer(p),as.double(Y), as.double(X),res=as.double(theta))
cat("\n");cat("time: ")
cat((proc.time()-t))
return(matrix(res1$res,nrow=M,ncol=p))
}
res<-ft()

dyn.unload("test/Rinterface.so")
